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We propose a soluble model to incorporate the nonlinear effects in the transition 
probabilities of the multiphonon Giant Dipole Resonances based on the SU(1,1) 
algebra. Analytical expressions for the multi-phonon transition probabilities are 
derived. For reasonably small magnitude of nonlinearity x ~ 0.1 — 0.2 enhance- 
ment factor for the Double Giant Resonance excitation probabilities and the cross 
sections reaches values 1.3 — 2 compatible with experimental data. 



1 Introduction 

Coulomb Excitation in collisions of rclatiidalkJons is one of the most promis- 
ing methods in modern nuclear physics Elfla'Q'El. One of the most interesting 
applications of this method to studies of nuclear structuEC is the possibility to 
observe and study the multi phonon Giant Resonances u. In particular, the 
double,Dipolc Giant Resonances (DGDR) have been observed in a number of 
nuclei QO'EI. The "hulk properties" of the one- and two-phonon GDR are now 
partly understoodid and they are in a reasonable agreement with the theoretical 
picture based on the concept of GDR-phonons as almost harmonic quantized 
vibrations. 

Despite that, there is a persisting diacrejpajicy between the theory and the 
data, observed in various experiments M'ElEl that still remains to be under- 
stood: the double GDR excitation cross sections are found enhanced bttrfectof 
1.3 — 2 with respect to the predictions of the harmonic phonon picturel 



This discrepancy, which almost disappears aLhjp'h brmilwdiug. energy, has at- 
tracted much attention in current literaturea,llj'tdlijli!]'tZl, E3£3lHj; among the 
approaches_jto resolve the problem are the higher order perturbation theory 
trpa trijprit P^l and studics of anharmonic/nonlinmiL aspects of GDR dynamics 
H,t3^L3p. Recently, the concept of hot phononslijEj within Brink- Axel mecha- 
nism was proposed that provides microscopic explanation of the effect. These 
seemingly orthogonal explanations deserve clarification which we try to supply 
here. 

The purpose of this work is to examine, within a a soluble model the role of 
the nonlinear effects on the transition amplitudes that connect the multiphonon 
states in a heavy-ion Coulomb excitation process. Most studies of anharmonic 
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corrections tSE^O concentrated on their effect in the spectrum E§E3. 

Within our model, the nonhnear effects are described by a single parame- 
ter, and the model contains the harmonic model as its limiting case when the 
nonlinearity goes to zero. We obtain analytical expressions for the probabili- 
ties of excitation of multiphonon states which substitute the Poisson formula 
of the harmonic phonon theory. For the reasonable values of the nonlinearity, 
the present model is able to describe the observed enhancement of the double 
GDR cross sections quoted above. 

Having in mind to show how analytical results follows from the nonlin- 
ear model, and to explain how the model works, we restrict ourselves here 
to its simplest version (transverse approximation, or SU(1,1) dynamics) and 
keep numerics up to minimum level. We postpone till further publications de- 
tailed numerical analysis and comparison with the data. Microscopic origins of 
the nonlinear effects (which are considered here phenomenologically) are also 
beyond the scope of our discussion. 



2 Excitation of multi-phonon GDR: Coupled-channels problem 

We work in a semiclassical approachS to the coupled-channels problem, i.e., 
the projectile-target relative motion is approximated by a classical trajectory 

ad the excitation of the Giant Resonances is treated quantum mechanically 
The use of this method is justified due to the small wavelenghts associated 
with the relative motion in relativistic heavy ion collisions. The separation 
coordinate is treated as a classical timejiependent variable, and the projectile 
motion is assumed to be a straight lineEI. 

The intrinsic dynamics of exciled. nucleus is governed by a time dependent 
quantum Hamiltonian (see Refs. Ei^j). The intrinsic state \il^(t) > of excited 
nucleus is the solution of the time dependent Schrodinger equation 



[Ho + V{t))] \i;{t) >, IV'(i) "^(*) 1^ > {-i^Nt) , 



. d\m > 

JV=0 

_ (1) 

where Hq is the intrinsic Hamiltonian and V is the channel-coupling inter- 
action. We let ?i = 1, c = 1. The standard coupled-channel problem is to 
find the expansion amplitudes aAr(t) in the wave packet as functions of 
impact parameter b where En is the energy of the state > specified by 
the numbers of excited GDR phonons A^. We assume the colliding nuclei to 
be in their ground states before the collision. The amplitudes obey the initial 
condition a^it — > —oo) — 5n.q and they tend to constant values as i — > ±oo 
(the interaction V{t) dies out at t — > ±oo). The excitation probability of an 
intrinsic state jA^ > in a collision with impact parameter b and the total cross 
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section for excitation of the state \N > (given by the integral over the imact 
parameter) are 



WNib) = \aNioo)\, C7N = 2n J bWN{b)db (2) 

(we use the grazing value bgr = 1.2{Ali^ + Alp^) as the lower limit). Here- 
after, the labels exc (sp) refer to the excited (spectator) partneiiin a colliding 
projectile-target pair. We neglect the here nuclear contributionEd to the exci- 
tation process. 

It is convenient to treat the coupled channel equations (Q) in terms of the 
unitary operator Uj (the interaction picture): 

i^Uiit) = Vimiit), Viit) = e^"°'V{t)e~'"'>\ Ui{t - -oo) - /, 
at 

(3) 

where the time-dependent Hamiltonian H{t) = Hq + V{t) that acts in the 
intrinsic multi-GDR states is given by Hq ~ LoNd, Nd = d'^dm and 

m 

V{t) = «i(t)[(£;i_i)t - (Sl+i)^] + vo{t){Ela)^ + Herm.Conj. (4) 

where Ell^ and El^ are the dimensionless operators actirig in the internal 
space of the multi-GDR states. The functions v are given inl3, e.g., 



Here, ruM and e are the proton mass and charge, Z , N and A denote the nuclear 
charge, the neutron number and the mass number of the colliding partners, 
7 = (1 — v'^)~^/'^ is relativistic factor, v is the velocity and the parameter p is 
the deal of the strength absorbed by the collective motion (usually assumed to 
be close to unity) EJ. 

In the harmonic approximation, the operators Ell^,Elm are given by 
the GDR phonon creation and destruction operators of corresponding angular 
momentum projection to, = This model of "ideal bosons" coupled 

linearly tp-the Coulomb field has well known exact nonperturbative solution 
(see, e.g. 113) for the excitation probabilities 

I harin \ 2 \ i harm 1 2 o I harm 1 2 , i harm 1 2 / a\ 

|a I = 2^ 1"™ I = 2|ai I + |ao | , (6) 

m=0,±l 
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i.e., the Poisson formula where the amphtudes ajjf are expressed in terms of 
the modified Bessel functions. At the coUiding energies sufficiently high, the 
longitudinal contribution (cx | is suppressed by a factor proportional 

to 7^^, see, e.g.,El. In the following, we will work in this "transverse approxi- 
mation" dropping the longitudinal term (the results are still qualitatively valid 
at lower energies). 

3 Nonlinear Model 

Our idea is to keep the spectrum of GDR system harmonic with the Hamilto- 
nian Hq = cuN. That is supported by the s^stematics of the observed DGDR 
energies, E2, which yields E2 — (1.75 — 2)a;lll. The conclusion on the weajL-aii- 
harmonicity in the spectrum follows also from theoretical considerationsE3p3. 

The transition operators El\ El that couple intrinsic motion to the Coulomb 
field can however include nonlinear effects: one can use the boson-expansion- 
like expression for them 




where the parameters Xi determines the strengths of the nonlinear contribu- 
tions. In particular, this could be a result of nonlinearities (anharmonicities) 
in the phonon Hamiltonian obtained in higher orders of perturbation theory. 
A reasonable way to take into account these nonlinear effects that we adopt in 
this work is to keep in the expansion (^ infinite series of type 

^^In ~ din + ^ '^m'^mi ^ '^m'^mi "^mi rfL2 + •■• = 

mi mim2 

= dl(l + 2xNdy^\ (8) 

where the single parameter a; > determines the strength of the nonlinear 
effects, i.e., the problem reduces to the harmonic oscillator with linear coupling 
when X — > 0. Nonzero values of x that we will consider here lead to a number 
of nonlinear effects that results in enhancement of the excitation cross sections 
for the double Giant resonances. 

Expression (^ , together with the standard coupled-channel formalism pre- 
sented in the previous section, composes our nonlinear model. Let us mention 
its advantageous points: 

(i) The model takes into account the nonlinear effects (not all of them by 
at least part of them) , nonlinearity is governed by a single parameter 

(ii) The model is exactly soluble beyond the perturbative theory 
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(iii) Its results reproduce correctly the magnitude of the enhancement fac- 
tor for the Double GDR cross sections and its dependence on bombarding 
energy 



4 Algebraic solution: SU(1,1) dynamics 

To solve the highly nonlinear problem (||) with and (||) we introduce the 
triad of operators D~^, D~ and that are given by 



D- - ^/k + -N,id+, - D+ = ^{dX, - dU)^2k 

(d+i - dU){d+i - + 2(2fc + Nd) 
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(9) 



where -D+ the Hermitean conjugate — (-D^)^ with Nd = d^id+i +c?l]^c?_i 
and k = (4a;) It is easy to check that they obey the commutation relations 
for the noncompact SU(1,1) algebra 



[D-,D"] - D-, [D+,D"] = -D+, [D-,D+] = 2D" 



(10) 



i.e., Eqs.(^) can be viewed a verision of the Holstein-Primakoff realization of 
SU(1,1) algebra. The parameter k = l/(4a;) is related to the Casimir operator 
of SU(1,1). 

The dynamics of the system, in the transverse approximation, can be 
expressed now in terms of the operators , D~ and (^ only. The 
interaction-picture evolution equatkm (^ and its formal exact solution, the 
time-ordered exponential (see, e.g.,Lj) than read 



Vk 



Vk 



Ui{t), 



Ui{t) = Texp I J dt'Vi{t') 1 (11) 



where the commutation relation [Nd,D^\ = ±D^ and the algebra ( |l0|) has 
been used in (|),(|) and (|). 

From purely mathematical viewpoint, the problem described by the last 
equation drops into the universality class of the systems with SU(1 



that can be analyzed by means of generalized coherent states for the 
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SU(Ll)j-algebra. (For other algebraic approaches to scattering problems, see 
Refsta,Eo|). 
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Due to closure of the commutation relations between the operators 
and that enter the exponential in Eg. (pi]), the time-ordered exponential 
can be represented in anothur equivalent form that involve ordinary operator 
exponentials only (see, e.g.,Ej): 



Uilt) — exp 



a{t) 



|a(t)P 



-it) 



D+ exp log {^1 - ' y' j - i(j){t) D° exp — '^^^ 

(12) 

and some time-dependent complex number a{t) (star means complex conju- 
gation) and real number (j){t) (phase) cB. The unknown functions a{t) and 
(j}(t) can be found from simple differential equations which relate them to the 
function vi{t) in the Hamiltonian H{t). These equations can be restored af- 
ter substituting the right hand side of Eq.(12) into the left hand side of the 
Schrodinger equation for the operator Ui{t) (11) and collecting the terms which 
have the same operator structure. Proceeding this way, we obtain, after some 
algebraic manipulations, from ( pT| ) with using the commutation relations (10) 
and applying the known formula formula for the operator exponentials 

1 



Y,X 



2! 



Y, 



Y,X 



the following Riccati-type equation for the complex amplitude a: 



vi{t)e"^' +Axvi{t)e~"^'a^ 
The phase (j){t) is given by a simple integral 



I — a 
dt 



(13) 



(j)(t) = (2/fc) / dtiRe{vi{ti)a{ti)e- 



In fact, (j){t) does not contribute to Wn and um- The simple nonlinear equation 
( p^ ) accounts for all orders of quantum perturbation theory for the problem 
Eqs.(p|),(P),(pl]). It is also seen from Eq.(p2[) that unitarity is automatically 
preserved within present formalism {U\ — Uf^). 

The expression for the amplitudes a n (t) which we are interesting in follows 
from (|l^) immediately after projection of the state 

\m) = ui{tm 

onto the intrinsic states with definite number of GDR phonons, A^. From 
Eq. (|l2|) , we have 



Ui{t)\0) =e-^'=^(*) (^1- i 



a(i)|' 



exp 



y/k 



|0). 
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From the last equation, we obtain the final expression for the amplitude |ajv(oo)| 
and the excitation probabilities for the excited states with N phonons: 



|ajv(cxD)| = (l — 4a;|a(a;)| 



Wn = |aAr(oo)p, 



1/2 



(4x|a(a;)|^ 



N/2 



(14) 



Here, the quantity a{x) is the asymptotic solution to the Riccati equation 
( |l3| ) at t — > oo subject to the initial condition a{—oo) = 0. The combination 
x^^|a(a;)| in (14) can be viewed as a "special function" of the two parameters. 



x^/'^F/uj and the adiabaticity parameter It can be easily tabulated by 
solving (13). The cross sections are then obtained from the usual formula (||) 
with using (p^. 

The harmonic limit of these results corresponds to the case x — > 0, when 
the nonlinearity disappears in the transition operators (^) and the coupling to 
electromagnetic field becomes linear. Then at a; ^ the last nonlinear term 
drops from the equation (|l^), and the amplitude is reduced to its harmonic 
value 



|a(x)| 



'±1 



that is given by _the modified Bessel function 
probabilities W 
harmonic results 



i^. The expression for the 
oo to the Poisson formula (||), thus the 



reduces at k 
are restored. 

At nonzero nonlinearity x > 0, the excitation probabilities Wn (^4|) for 
multiple GDR {N > 1) turn out to be enhanced as compared to their values 
in the harmonic limit W^""""", as illustrated in Fig.l. 

The deviation in the excitation probabilities for the A^-phonon states Wn 
( p!^ ) from their harmonic values W^""™ (^) (the enhancement factor) is given 
by the ratio 



Wn r(2^ + 7V) (l-4x|a(x)|2)^ |a(x)|2^ 



TJ/ Alarm 

y'' N 



,-2|a5"'"' 



Jiarm\2N 



> 1. 



(15) 



The following ratio is especially representative (see (|6|),(p^)) for the two- 
phonon case 



i?2 = 



W2 I Wi 



yynarm 



I wl' 



= (l + ^'^) uharm|2 - ^^^^ 
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Figure 1: Effect of the nonlinearity on the Af-phonon excitation probabilities (schematic 
plot). Nonlinear model values (dark circles connected by the solid line) as compared to 
the harmonic model (Poisson) results Wharm (empty circles connected by the dashed line). 



The first factor in this expression, 1 + 2x, results from the kinematic enhance- 
ment of the transition probabilities due to nonlinear effects considered here. 
The last factor in ([l5|),(p^ results from dynamical effects caused by linearity 
which are incorporated in the asymptotic solution of the nonlinear equation 
(p^). Unlike the first factor that depends on x only, the second one depends 
on the bombarding energy and it gives rise to additional enhancement in low 
energy domain. The same is valid for the emnhancement factor in the cross 
sections, r2, that is given by the formula 

oo 

2tt J bW2db 

harm oo ■ U'J 

2 27r / foVF2'"'™d6 
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5 Results and discussion 

The interesting feature of these results is that the enhancement factor (|l^) , (^6|) 
is much more sensitive to the bombarding energy than to the parameters of 
the spectator partner. This is just what has been observed in experiments: 
the values of r^c^^ found foij-.DGDR in ^^^Pb projectile using different targets 
^^^Sn, ^^^Ho, ^^^Pb, "^^^Ua are very close to each other and they correlate, 
within the error bars, with the value 

r2f"*P6) ~ 1.33 , high energy, 7 ~ 1.7 (18) 

(bombarding energy e ~ 640 Mev/per nucleon). 

The same picture was found in experiments on Coulomb deaintegration 
of ^^'^Au target using various projectiles "^^Ne,^^ Kr,^^'^ Au,'^^^ BiB. Also, the 
similar conclusion of nearly constant value of r2 in ^^^Pb target with various 
projectiles has been made in worktSI but for low bombarding energy e ~ 60—100 
Mev/per nucleon. In this case, 

r2(^°^P6)~2 , low energy, 7 -1.06 -1.10. (19) 

According to the results of calculation within present nonlinear model, 
these enhancement factors would correspond to the nonlinearity parameter x 
equal to 

x{^°^Pb) ~ 0.16-0.20 

i.e., to a reasonably small nonlinearity. 

Below, we present the exact results for the cross sections calculated ac- 
cording to Eqs.(|^), and ([l^) , (|l7|) and with solving Eq.(|T^) numerically. 

The dependence of the enhancement factor r2 — uila^""™ (|l7| ) for the 
DGDR excitation on the strength of the nonlinearity x is shown on Fig. 1. 
for the process ^'^^Pb +2*'^ Pb (we use the case of bombarding energy e — 
0.64G'ey/per nucleon). It is seen that the enhancement factor drops to unity 
at big values of k (harmonic limit) and grows at stronger nonlinearity. The 
scaling value of r2 is also shown for comparison. 

In view of the difference between the enhancement factor values at high 
and low energies (Eqs.( |l8| , |l9|) ), it is interesting to trace energy dependence of 
the enhancement factor. 

Deviations in r2 from the straight line 1 + 2a: occur at both low and high 
energies. At 7 — s- 1, adiabatic approximation is valid, and this yields to \a\ > 
|Q,ft,arm|_ q;^}jys^ i?2 > 1 + 2x and 7'2 > 1 + 2a;. By contrast, at higher energies, 
the dynamical nonlinear effects tend to reduce the magnitude of thus 
|a|/|Q;5""'™| < 1, and i?2 < 1 + 2x. To sum up, the enhancement factor 
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Figure 2: Enhancement factor r2 = 0-2/0-2'"''" for the Double GDR excitation in 208p^^208 
Pb process at bombarding energy e = 640 Me V /per nucleon as a function of the inverse 
nonlinearity parameter k = (symbols, solid curve is to guide the eye). Dashed curve 
corresponds to the scaling r2 ~ (1 + 2x) = ^'^^ ■ 



for the DGDR excitation cross section, r2 = a2/cT2°'^™' drops from 2 — 2.5 
(for low bombarding energies e ~ lOOMeF per nucleon) to 1.2 — 1.3 (for 
e 640 — IQQMeV per nucleon) -while fixed value of nonlinearity x is used. 

On Fig. 2., -we plotted the value of the enhancement factor calculated nu- 
merically for the case of "^^^Pb +208 P5 process. The magnitude of nonlinearity 
is kept fixed x = 0.19. 

It is seen that reasonably small nonlinearity reproduces correctly the ob- 
servable value of the enhancement factor and its energy dependence. These 
results are is in correspondence "with the experimental]^ observed trends and 
-with microscopic models based on Axel-Brink conceptllj. 

To conclude, -we presented here a simple model that accounts for the non- 
linear effects in the transition probabilities for the excitation of multi-phonon 
Giant Dipolc Resonances in Coulomb excitation via relativistic heavy ion col- 
lisions. The model is based on the group theoretical properties of the boson 
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Figure 3: Enhancement factor r2 = cr2 / cr!^'"'"^ for the Double GDR excitation in 208p{|^208 
Pb process as a function of rclativistic factor 7 (symbols, solid curve is to guide the eye). 
The value of the nonlinear parameter x is kept to be equal to a; = 0.19. The scaling value 
(constant 1 + 2x = ) is shown by dashed curve. 

operators. It allows to construct the solution for the dynamics of the multi- 
phonon excitation within coupled-channel approach in terms of the generalized 
coherent states of the corresponding algebras. The exactly solvable harmonic 
phonon model appears to be a limiting case of the present model when the 
nonlinearity parameter x goes to zero. The main advantages of the limiting 
harmonic case (unrestricted multiphonon basis, preservation of unitarity and 
analytical results in nonperturbative domain) remain present in our nonlinear 
scheme. Therefore, the model can be viewed as a natural extension of the 
harmonic phonon model to include the nonlinear effects in a consistent way 
while keeping the model solvable. 

The probabilities to excite double-phonon GDR appear to be enhanced by 
a factor (1-|-2x)(|q:(x)P/| a^'"'"'" | ^ ; that results in the cross section enhancement 
factor. This can be viewed as a hint that the discrepancy between the measured 
cross-sections of double GDR and the harmonic phonon calculations can be 
resolved within present nonlinear model by means of using an appropriate 
value of the nonlinear parameter x for a given nucleus. The experimental 
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values of enhancement of with respect to the harmonic results for the excited 
208p^ nucleus are almost insensitive to the details of the collision process. The 
enhancement factor drops as the bombarding energy grows. This is consistent 
with the data and gives results similar to those recently obtained in a possibly 



different context, with a theory 
and the Brink- Axel mechanism 
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the concept of fluctuations (damping) 
. It would be certainly worthwhile to 



establish possible connections between the two approaches. 
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